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Abstract 

We give a self-contained exposition of the combinatorial solution of quantum mechanical 
systems of coupled spins on a one-dimensional lattice. Using Trotter formula, we write the 
partition function as a generating function of a spanning subgraph of a two-dimensional lattice 
and solve the combinatorial problem by the method of Pfaffians provided the weights satisfy 
the so-called free fermion condition. The free energy / = J, h) and the ground state 
energy cq = eo(J,h) as a function of the inverse temperature /3, couplings J and magnetic 
fields h, for the XY model in a transverse field with period p = 1 and 2, is then obtained. 



1 Introduction 

A quantum mechanical system of coupled spins on a one- dimensional lattice A G Z with |A| = 
sites is described by Hamiltonians of the form 

HN{a) = -Y,hi,i+i (1.1) 



where 



h,, = Jl^ at + 31,0^0] + Jl^ at + ^{at + a]) . (1.2) 

For each z G A, af,a = x,y,z, are Pauli matrices and periodic boundary condition is imposed, 
o't+N = ^i- The family of couplings J = {(Jj^j+i, Jli+i, Jli+i)}i£A and magnetic fields h = {/ijjjgA 
we shall considered here are assumed to be periodic of period p and, for simplicity, only p = 1, 
Jfi+i = hi = h and p = 2 arrays, Jf^+i = Ji, hi = hi ii i is odd and Jfi+i = J21 hi = /12 if « 
is even, will be treated. The 2;-component couplings Jti+i will also be set iqual to zero at certain 
moment. 
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One-dimensional quantum systems have been extensively studied, as compared to the higher 
dimensional counterparts, because some of their physical quantities can be explicitly calculated. In 
the present work we consider the free energy function, / = /(/?, J, h), and the ground state energy, 
Co = eo(J,h) = lim /, where 

/=lim^lnZ^. (1.3) 

The partition function 

Zn = Tr e-^^^ = Tr p^v , (1.4) 

is given by the trace over a density matrix pat, defined in the 2^-dimensional vector space (3)igA C^, 
and (3 is the inverse temperature. 

Among the available solution methods, we distinguish two types depending on whether one 
approaches the partition function Zi^ directly or by means of the eigenvalue problem of H^- 

The eigenvalues are more easily accessible when, by a canonical transformation, Hj^ can be 
written as a free fermion system. If J and h have period one with = 0, then (jl.lj) is unitarily 
equivalent to [K] 



n 



N 



-t J2 {4 Q+i - c, + r (cl 4, - c, c+i) +h(i-2cI Ci) } , (1.5) 



where c| and q are the creation and annihilation operators for electrons at the site i G A,t = J^+J^, 



r = (J^ — J'^)/t, H = h/t and the boundary condition terms has been omitted. Lieb et al |LSMj 
and Katsura |Kaj have calculated the spectrum of the XY-mode\ ( ()1.5j) with H = 0) and XY- 
model in a transverse field (jl.Sp . respectively, obtaining the free energy function and the ground 
state. Subsequently, the XF-model with alternating interactions and magnetic moments has been 
solved by Perk et al fPHZ] . 

The partition function is said to be approached directly if it is mapped into a problem of counting 
spanning subgraphs in a lattice. Direct approach does not require Hj^ to be diagonalizable but in 
order to solve the underlying problem certain conditions have to be met. The basic ingredient for 
this approach is the Trotter formula [Tj, 



eV^^+^^)= lim fe^V-e^^/-r , (1.6) 

m— >oo 



valid for any finite-dimensional matrices K^, (see e.g. jS] for extensions and applications). The 
partition function ()1.4|1 can thus be mapped into a classical two dimensional lattice spin system as 
explained in detail in the following sections. 

There is another relationship between chains of quantum spins and the transfer matrix of two- 
dimensional classical systems which has been discussed in the excellent review by Kasteleyn |Klj . 

The partition function has been directly approached by Suzuki and Inoue |SuIj who, using 
Trotter formula, have mapped the period-one Xy-model in a transverse field into a period-two 
8-vertex model with two sets of weights, {wi}f^i and {w^}f^l, satisfying the free-fermion condition 

WiW2 + W3W4 = W5W(i+W7Ws. (1.7) 

Equation ()1.7p is the very fact that allowed Fan and Wu |FWj to solve the 8-vertex model by the 
Pfaffian method (see jHLWj for the period 2 model). 
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The purpose of this paper is to give a step-by-step account of the combinatorial method of solv- 
ing chains of interacting quantum spins. In particular, Katsura's results are shown to be recovered 
by mapping the XF-model into a period-one 8-vertex model with weight function satisfying (|1.7|) . 
For this we use the one-to-one mapping established by Barma and Shastry |BSj in the context of 
XYZ-model. Due to the careful treatment of our exposition the method is extended to any period 
p. To our knowledge, we calculate for the first time the free energy of a period-two XY-mode\ in 
a transverse magnetic field by the method of Pfaffian. 

This paper is part of a program aiming at alternative solutions to the ground state energy of 
Heisenberg fll.ip models. The method based in the spectral analysis of Hn is amazingly far more 
developed than the direct approach. The ground state energy of Hubbard (LWj and Heisenberg 
models |Y Yj can be obtained by an Ansatz to the eigenfunctions of if a? called Bethe-Ansatz. The 
need for alternative solutions is more demanding since the "string hypothesis" (completeness of the 
Bethe-Ansatz eigenstates) remains unresolved. 

Last but not least, the interest on the Pfaffian method has been renewed in view of a recent 
result by Pinson and Spencer |PSj on the universality of critical d = 2 Ising model. 

The layout of this paper is as follows. In Section |21the Hamiltonian of the quantum spin system 
(jl.lj) is mapped into a two-dimensional classical spin system. This in turn is mapped, in Section 

01 into a spanning subgraph problem which is solved in Section |3] by the Pfaffian method. The 
evaluation of the determinant is presented in Section 14.11 and free energy functions for periodic 
systems with periods p = 1 and 2 are determined in Sections IHl and [71 respectively. 

2 Trotter Formula 

In this section the Trotter formula will be used to write the partition function (jl.4|) as the partition 
function of a two-dimensional lattice x Z"* of a classical spin system in the limit m oo. 

Trotter's product formula ()l.fi|l is a mathematical tool which has elucidated the "path integra- 
tion" of a quantum mechanical particle in a potential (see [S] and references therein for an account- 
able survey). The application of this formula for coupled quantum spin systems was originally 
carried out by Suzuki ^Su] but we shall here follow a slightly modification due to Barma-Shastry 
|BSj . Before working it out, let us define the Hamiltonian system precisely. 

To each site i G A it is associated a Hilbert space Hi = C^. The vectors e+ = (g) and e_ = Cj*) 
form a base of Hi and represent the two possible states of a spin Sj: Sj = -|- if the spin at i is up 
and Si = — if it is down. 

The vector space H = Ha where the Hamiltonian operator acts is a tensor product space, 

H = ^Hi. (2.1) 
Under the above interpretation, each of the 2^ vectors of the form 

represents a possible spin configuration s = {sijjgA of the whole system. Note that the collection 
B = Ba = {cs} of vectors Cg := e<j^ ® e<j2 ® • • • ® e^j^ with s G {+, —}^ forms an orthonormal base 
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of Ti. So, if Es = Cs (cs)^ denotes the (orthogonal) projector of TC into the Cg direction, we have 

Yl ^^=^' (2.2) 

where J=/(8)/(g>---(8>/is the identity matrix: Jf = f , for all G 7^. 
A base of operators in can be chosen to be the Pauli matrices 

^'-{"1). ^-CoO- ^-(J-O- 

and the identity / = ( ^ ? ) • For later reference, we shall also consider the "raising" and 



'lowering" matrices 
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whose operations over the base vectors are e_ = e+, e+ = e_, and e+ = cx^ e_ = 0. 
The spin operator erf in (|1.2p acts on the Hilbert space Ti and has the form 

(Tf = /®---®/®(T"®/®---®/ (2.5) 

with 0"° acting on TYj (i.e. it is located at z-th position). Here a runs over x, y, z, 0, + and — with 
= I. The Hamiltonian is composed of pair interactions between nearest neighbor sites, 

< af+i = / (g) ■ ■ ■ O / ® a" (g) (g) J (g) ■ ■ ■ (g) / , 

which acts as the identity operator in Tik for all k except at the positions k = i,i + 1. The 
type of interactions considered in ()1.2|) have pairs {a, (3) equal to (x, x), {y,y) and {z,z) which will 
subsequently be written in terms of (+, +), (+, — ), (— , +) and (— , — ) interactions. 

Note that two consecutive terms of the Hamiltonian (jl.lj) do not commute, 7^ 0, 

due to the fact that 



does not vanish if /? 7^ 7. This expresses the quantum mechanical nature of our spin system. We 
shall now explain how Trotter's formula resolves the non-commutativity character by converting 
quantum spin variables } into infinite many copies Si, S2, . . . of "classical" spin variables Sj = 

For N = 2n let us decompose the set A = {1,2, . . .2n} = A°'^'^ U A''^^'^, according to whether 
the site i is odd or even, and write 



H2n = H:'"' + Hr\ (2.6) 

where }{°'^^^^^^^^ jg given by with A replaced by A"'^*^^'^^'''^^. Because we have eliminated con- 
secutive pairings, the terms which compose commute with each other. Note, however, 
that [Hf'^, i/^™"^] 0. 
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Using Trotter's formula p.6|) we have 

P2„ = lim (e-^''"-''/"^ e-^^"^^"/-)™ ^ lim {p^^pTr^r ■ (2-7) 

m— ►oo V / m— >oo ' ' 

The partition function ()1.4|1 can thus be written as 

Z2n = hm Zn,m (2.8) 
m— »oo 

where 

7 — Tr n°'^'^ even odd even |'9 Q^l 
^n,m -LJ- Hn,m h'n,m l-'n,mHn,m ■ 
V ' 

m— times 

We are going to make use of our base B by inserting the spectrum decomposition of the identity 
fj2.2p between all pairs of density matrices pO<id j^even ^^^^ ^even j^odd_ other words, we want to 
represent the density matrix p according to the base B. To make this representation clearer we 
shall first consider the special case Zi i. 

Dropping, for simplicity, the subindex of p, we have 

= Tr p°'^^ p'''"'-^ = Tr p"'^'^ Ip''"^'^ ^=Y1 -^2' P-'T (2-10) 

s,s' 

where ps,s' = {^sY P'^s'- For n = 1, A = {1, 2} has only two sites and 7f is a four dimensional space, 

s e {++, H — , — h, }, leading to a 4 x 4 density matrix p. In view of the periodic boundary 

condition, p'^^'^ differs from p^"^^"^ by the order of the spins at the sites {1,2}. We shall keep it in 
mind and drop their super-index too. As we shall see, the alternating position of the odd-even 
density matrix in the product ()2.9p will lead to a chess-board pattern. 

To compute the 4x4 density matrix p, the algebraic properties of the matrices ()2.Hj) and ()2.4|1 
have to be used. Let us first write /ii,2, with = 0, in terms of cx'^,a~ and cr^. From p.2|l and 
flTl) . we have 

/ii,2 = (r - jy) {at at + ar a,) + (r + jy) {at a, + ar at) + ih/2) {a^ + a^,) , (2.11) 
where Jf 2 — J21 — J'^ ^"^^ hi = h2 = h. 



Writing rj = {J^ — J^), k = {J^ + J^) and C, = ^/rf + h'^, the following relations 



(M''e. = <' 'j^,:--' :_ ' ' (2.12) 
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(M''"'e. = <( ViT-^' ' (2.13) 

hold for any integer number k>l. 
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Taylor-expanding p = e^^^'"^ and using equations p.l2|l and p.l3|l . combined with the orthog- 
onal relation (es)"^es' = 5s,s', gives 



/ cosh/3C+ (VC)sinh/3C (VC) sinh \ 

cosh /3k sinh/^K 

sinh/3K cosh /5k 

\ (r//C) sinh /?( cosh /?( - {h/Q sinh PC J 



(2.14) 



Note that p is a symmetric matrix and Zi^i = Tr p^. 

Now, we want to include the term into the above calculation. The simplest way to deal with 
this case is to introduce a slightly modification of the Trotter formula |KMj . 



K0+K1+K2 



lim (< 



,Ko/in Ki/m K2/in 



)' 



As in |BSj we separate the J^-dependent terms of and divide the remaining terms according 
to the even-odd partition of A, 



H2n = m^ + Hf'' + Hr\ (2.15) 

where jg f|2.fjj) . We repeat the steps ()2.7j) - ()2.9|) with the modified Trotter formula. 

For the special case Zi 1, we have 

= Tr p° p"'''' p^^^'^ = Tr (p^) p"'^^ p^™"^ (p°) = ^ p° J p^™" , (2. 16) 



where p°'^'^ = (p°)^''^ P = P with p given by ()2.14p . 



( e2/5^^ 



V 








(p)"^ and 








\ 






Note that p = (p)^. 



(2.17) 



To extend the calculation of 1 to ^„,m with n,m > 1, we may think of a rectangular 2m x 2n 
lattice as a "chess-board" with n ■ m black squares. To each of these black pieces labeled by p (of 
plaquetes) we associate a matrix pp given by (|2.14p with jS replaced by /?/m; p is a transfer matrix 
between the state of spins at the up corners of the black square p and the state of the spins at the 
bottom corners of p. If the up and bottom corners' labels of p are respectively v, w and v', w', we 
denote the matrix elements of pp by p(st,, s^; s„/, s^.') and equation ()2.9|) represented in the base B 
can be written as 



m In 



Zn,m — W I \\_P{s2i-l,2j-l^ S2ifij-l] S2i-lfij, S2i^2j) P{s2i,2j^ S2i+l,2j'-, S2i,2j+l^ S2i+l^2j+l) 



(2.18) 
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with periodic conditions assumed in both directions. For ()2.18|) . note that the matrix elements of 
podd(cven) between two successive spin configurations, = {sk^i}lZi and S;+i = {sk,i+i}l^i, factor 
out 

n 

{Pn%) s2j_^^S2j ^ nP('52i-l,2i-l, S2i,2j-i; S2i-l,2j, ■52j,2j) , 
i=l 

and analogously for p"^^^, due to the even-odd decomposition. 

As a consequence, two matrices pp and pp' are coupled if p and p' share the same vertex v 
[pCip' = v). Each matrix labeled by a plaquete p has a pairwise interaction with four other matrices 
labeled by the nearest-neighbor black plaquetes of p. We shall use this fact in the following section 
to write ()2.18|) as a counting problem in a n x m rectangular lattice (wrapped into a torus and 
rotated 45°) whose vertices are centers of the black plaquetes of the original lattice. 



3 The Vertex Model 

In the previous section Trotter's formula has been used to write the partition function of a 
one-dimensional quantum spin system in = {!,••• , 2n} as the m — > cxo limit of a sequence 
{Zn,m}m>i of partition functions of two-dimensional classical spin systems in An,m = ^21% x ^2-01 
with four spin interactions mediated by the matrix p (see equations ()2.8j) . ()2.9j) and ()2.18|l ). The 
mapping we shall now discuss is due to Barma and Shastry |BSj . 

We consider A„ „j wrapped in a torus as a "chess-board" and define a n x m regular lattice A* ^ 
with vertices at the center of the black plaquetes of A„ „ and edges connecting nearest neighbor 
vertices (see Figure 1). Note that the edges of A*^ are 45° with respect to edges of A„ ,„- The 
remarkable property in this construction is the one-to-one correspondence between vertices of A^^m 
and edges of A* ^. 



Figure 1: The regular A„ „i and dual A* ^ lattices 
Let us define the "dual" relation * : A„_m — ^ A* defined by 

{hiT = e, 

if the the edge e G A* ^ crosses the site with coordinates (i, j) G A„^m, and 



if the vertex v G A* ^ is the center of a black plaquete p e A 
a correspondence between classical spin configurations {sij 
the two-state edge configurations in A*_^. We say that the edge e 



n,m- 

i = 1 



This allows us to establish 
. . . , 2n; j = 1, . . . , 2m} and 
is occupied or vacant 



according to whether the value of Sij is 1 or —1. 

Now, the density matrix p at the plaquete p, whose values depend on the spin configuration of 
its four corners, is mapped, by the "dual" transformation, into a vertex function w ai v — p* whose 
values depend on the configuration of the four incident edges. In the period 1 system there is 
only one vertex function w for all v e A*^^. This leads us to the following definitions. 

A vertex configuration is an assignment ^ : A* _^ — > {1, . . . , 8} which associates to each vertex 
V one of the 4-incident-edges configurations shown in Figure 2. There are actually 2^ = 16 
possible assignments of occupied-vacant incident edges but only 8 of them (shown in Figure 2) are 
relevant to the problem at hand. 



X 



5 

X. 
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X 



Figure 2: Possible assignements of 

A vertex configuration ^ is said to be compatible if for each edge e with end points v and v', 
and $,y' have the same assignment for the common edge e. From here on, by vertex configuration 
^ we always means a compatible configuration. It is worth noting that a vertex configuration ^ 
corresponds to a covering of A* by closed polygons (see Figure 3) . 
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Figure 3: A vertex configuration ^ 
The weight of a vertex configuration ^ is given by 



(3.1) 
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where, in view of fj2.14|) . ()2.16|) . the "dual" mapping and the label chosen in Figure 2, 

m Q m 



W3 = W4 



m ( m 
g-/3JV-ginh^, (3-2) 

W5 = ujq = e"^"^^''™ cosh — , 

W7 = Ws = e"'^"^^/™^ sinh — , 

C m 

where (, rj and k are given following equation ()2.11|) . Note that only the eight vertex configurations 
shown in Figure 2 have non-vanishing vertex function. 

With these definitions, we conclude 

Proposition 3.1 The partition function ^2.1^) can be written as 

Zn,m = ^w{0 , (3.3) 

where the sum runs over all compatible vertex configurations ^ with weights given by i3. 1]) and 

At this point, one can verify whether ()3.2j) satisfies the free fermion condition stated in ()1.7j) . 
The equation 

w^W2 + wsw, - w,w, - w-jw, = {e'^''/^ - e-^^'^l^) (^1 + ^ sinh^ 
vanishes if and only if = 0. This will be assumed in the next sections. 



4 Pfaffian Solution 

In this section the partition function Z„ m, written as the generating function of a spanning subgraph 
of A* = A*^ in ()3.3|) . will be solved by the Pfaffian technique. For this, we introduce some 
terminology. 

Let ^ be a graph which we may think to be a "decoration" of A*. A dimer covering D oiQ ys, 
a spanning subgraph consisting of non-overlapping close-packed configurations of elements formed 
by an edge and its two terminal vertices (dimer). Note that no vertex is left vacant in D. 

In the early 60's the problem of enumerating all dimer coverings has been solved for certain 
lattice graphs (see (K] for a review). Right after this. Green and Hurst [GHj showed that the 
counting dimer coverings can be useful to solve the partition function of two-dimensional Ising 
models (see also |GHj ). Since the Ising model has to be written as the generating function of a 
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spanning subgraph of the form ()3.3|) . the method we now start to explain can be apphed directly 
to this equation. 

Given a collection of positive numbers z = {ze}e£g indexed by the edges of Q, let r(z) be the 
generating function for the dimer covering problem on Q, 

T{z) = J24D), (4.1) 

D 

where 

ziD) = n 

eG-D 

is the weight of a dimer covering D with the product running over edges covered by the dimers of 
D. 

If {1,2, . . . , k} is an arbitrary enumeration of the vertices of Q let A = {a^y} be a skew- 
symmetric matrix of order k whose elements a^y = —ciyx satisfy 

{Ze, if e = (xy) e Q , 
(4.2) 
, otherwise . 

The sign of a^y is chosen according to a given orientation of the edges of Q: a^y > if e is oriented 
from X to y. 

Our ability to calculate the partition function ()3.3p depends on whether the graph Q 

obtained by decorating A*, and the associate weight function, agree with the hypotheses of the 
following results. 



Proposition 4.1 Let A be a k x k matrix defined as above with k even. Then 

detA = (PfA)2, (4.3) 
where the Pfaffian, Pf A, of A is defined by 

PfA = 5^£p^a,,, 



with the sum over all partitions of {1, 2, . . . ,k} into pairs; is a permutation ji j2 ■ ■ ■ jk such that 
lii i2b3 isl ■ ■ ■ \jk-ijk\ is a description of the partition tc and ep^ is the signature of P.„ 



For a proof of Proposition 14. II see jKj. To understand, however, equation ()4.3p and get a feeling 
of the next statement, let us observe that if we superpose a dimer covering D on another dimer 
covering D' we obtain a circuit covering of Q. For, starting from any vertex of Q we walk along 
the edge in D till we reach other extremity; then, walk along the edge of D' and along the edge 
of D after the extremity has been reached and so on. Continuing this process we eventually arrive 
at the starting vertex when we have to choose a different one and start again this process till all 
vertices of Q have been covered. It is a matter of fact that det A is the generating function of circuit 
covering of Q. Note that 

det A = ^ Ep ai^p(i) ■ • ■ ak^p(k) , 
p 
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and all permutation P = {-P(l), • • • , -P(fc)} of {1, ... , k} can be decomposed into cycles (Ci, . . . , Cg). 
Since by the superposition of D on D' with D, D' running over all dimer covering one gets all circuit 
covering of Q, it is tempting to identify detA with T'^iz). As we shall see, this requires to choose 
the orientation of Q properly. 

Proposition 4.2 // the number of contra-oriented steps, with respect to the chosen orientation of 
Q , performed along each circuit generated from the superposition of D on D' is odd, for any two 
dimer coverings D and D' of Q, then 

r(z) = |PfA|. (4.4) 

By imposing that the parity of all circuits to be odd det A becomes a sum of positive terms^ 
and consequently the sign of every term in Pf A remains the same implying Proposition 14.21 (see 
ref. jK], pg. 89, for a proof). We shall call an orientation of Q with such property admissible. 

Combining the two propositions we arrive at 

r(z) = VdetA. (4.5) 

It thus remains only to write ()3.3|) as a generating function of dimer covering of a properly oriented 
graph Q and evaluate the determinant of the associated skew-symmetric matrix A. 

We begin by choosing the appropriated "decorated" graph Q. We take Q as being the terminal 
graph (A*)-^ (see Figure 4) obtained by the replacement of all vertices f G A* by a cluster of 4 
vertices all connected with each other by single edges (a complete graph K4) denoted by D^. Each 
vertex of Q is an extreme or terminal of an incident edge e G A* a property by which Q is named. 
The edge e G ^ will be called internal if it belongs to n^, for some v and external otherwise. 
The external edges are identified with the edges of A*. Note that the graph Q is not planar due to 
the crossing of diagonal edges of n^, which prevents, in general, that an admissible orientation can 
be found leading the Pfaffian method useless to compute r(z). As we shall see, it is exactly this 
problem that will help us to solve another one. 



Figure 4: Construction of the terminal graph (A*)^ 

What makes a terminal graph significant for our problem is that we can associate to each 
(compatible) vertex configuration ^ in A* a dimer covering D in the terminal lattice Q as follows: 

-^Note that the signature ep of a permutation P is given by (—1)* where s is the number of cycles (Ci, . . . , Cs) 
in P. If the parity of all circuits is odd, we get (—1) from the product of matrix elements along each cycle Ci 
compensating the signature. 
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an external edge e G -D if the edge e G A* is occupied in the configuration ^. A hst of vertex 
configurations at vertex v and its associated dimer covering at is shown in Figure 5. The 
problem here is the degeneracy of dimer covering associated with ^„ = 2. 




5 6 7 8 



Figure 5: Vertex and their associated dimer configurations 

It turns out that there exist a (non-admissible) orientation of Q such that, if a dimer configu- 
ration D' differ from D by the substitution of a crossing dimer covering at by a non-crossing 
one, the superposition of D' on D yields a circuit in with even parity and a minus sign is given 
to this contribution. In addition, this orientation is such that all circuits, excluding the ones with 
crossings and those winding around the torus, generated by the superposition of D on D', with D 
and D' two dimer coverings in Q, have odd parity. The orientation of Q with these properties is 
shown in Figure 6 for a given 




Figure 6: The orientation of a cell 

As demonstrated in ref. [K], pages 98-99, it is the crossing dimer configuration which gives rise 
to the minus sign and Proposition 14.21 holds if z{D) is replaced by z{D) = {—ly z{D) in equation 
M.lj) where c is the number of v's such that n^, is cross-covered in D and free boundary condition 
is imposed in Q. 

As An,m is wrapped on a torus. A* ^ may have many different topologies depending on the 
values of (n, m). The lattice A*^ will be wrapped on a torus only if m is a multiple of n. We 
shall in the following assume that m = kn+ 1 for some integer k which makes A* ^ to be helically 
wound on a torus and avoids circuits winding around Q = (A* „j)"^ (see McCoy and Wu |McWj for 
the dimer covering problem on a torus). 
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Once the orientation of Q has been resolved, we are in position to write ()3.3|) as a generating 
function r(z) of dimer covering of ^ = (A*)-^. We observe that the degeneracy in the association 
of ^t, = 2 resuhs from the dimer covering of using only internal edges. For all dimer covering 
associated with C,v 7^ the association is unique and there are at least two external edges covered. 

Now, let Q denote the class of equivalence defined as follows. Two dimer covering of Q, D and 
D\ are said to be equivalent, D ~ D' , if they give rise a configuration ^ under the association made 
in Figure 5. Then equation ()4.H) can be written as 

r(z) = $^«(0, (4.6) 

where 

u{0 = <D) . (4.7) 

As Fan and Wu |FWj (see also |HLWj ). if in addition we set Ze = 1 for all external edges e G ^, 
then u{S,) factors out, 

^(0 = n (4-8) 

where Mg„ depends only on {zejeen^ and is explicitly given by (using the label of Figure 5) 



1, 


if 


e = 


= 1, 


2^(14) 2(23) + 2(34) 2(12) - Z(i^)Z(^24) , 


if 


e = 


= 2, 


2(24) , 


if 


e = 


= 3, 


2(13) , 


if 


e = 


= 4, 


2(34) , 


if 


e = 


= 5, 


2(12) , 


if 


e = 


= 6, 


2(23) , 


if 


e = 


= r. 


2(14) , 


if 


e = 


= 8. 



(4.9) 



In order to compare equation ()3.3p with ()4.6|) . we identify 'u(^) with the normalized vertex 
function w{^)/wi: u^^ = w^^/wi. Equation ()3.1|) then reads 

^n,,n = <'"5Z«(0=<™-r(z). 

Note that z = {z^} has to be specified only for the six internal edges of for a given v (recall 2e = 1 
for all external edges). Moreover, note that the last six equations of ()4.9|) determine completely the 
correspondence between the vertex function u{C,) = w{^)/wi and the weight z since, by substituting 
those equations into the second equation, we have 

U2 = + urUs ~ U3U4 = K = — , (4.10) 

wf Wi 

provided Wi, . . . ,Ws satisfy the free-fermion condition (|1.7|) . 

In view of fl4.5|) . the conclusion of this section may be summarize by the following result. 
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Proposition 4.3 Let A be a [Anm x 4nm) skew-symmetric matrix defined as in with the 

sign determined by the orientation of Q as in Figure 6. If Ze = 1 for all external edges and the 
equations = w^/wi, ^ G {3, 4, ... , 8}, given by j^.^j ) define for all internal edges e E Q, then 

Z„,^ = <"^v/d^. (4.11) 



5 Evaluation of the Determinant 



In this section we evaluate ttie determinant of tlie Anm x Anm matrix A, with entries given by 
fl4.2p . as described in Proposition 14.31 

We begin by block-diagonalizing A by the Fourier method. Let us first remind some facts about 
Kronecker (or tensor) product. If M = {rriij} and = {riij} are kxk and / x I matrices, respectively, 
the Kronecker product of M and N, M ® N, is a. kl x kl matrix given hj M ® N = {rriijN}. If 
Ml ® Ni and M2 ® are two matrices as above, their product is given by 

(Ml (g) A^i) (M2 ® N2) = Ml M2 A^i N2 . (5.1) 

The Kronecker products do not commute, M ®N ^ N ® M, but there exists a permutation matrix 
P such that P^{M^N)P = N(g)M. So 

detM ®N = detN ® M = det (/ ® M) det {I ^ N) = (det M)' (det N)'' . 



We want to write A as a nm x nm matrix whose elements are 4x4 matrices, called inner 
matrices. Since the entry axy of A is nonvanishing if (xy) = e is either an edge at or connecting 
nearest neighbor ones, there are only three relevant inner matrices. The inner matrices are labeled 
by the enumeration of vertices of □„ and the sign of their entries determined according to the 
orientation of Q (see Figure 6). 

Let U be the skew-symmetric "adjacency matrix" of □^, defined by the last six equations of 

(SD, 



u 



and let X and Y be defined by the external edges bridging 
southwest directions, respectively. 










—U4 


-Us 


\ 









U7 


Us 






U4 


-uj 





-M5 




v 


Us 


-U3 


M5 





/ 



(5.2) 



and in the southeast and 



X 



/ \ 



10 

\ / 



and 



Y 



/ \ 

1 



\ y 



(5.3) 



The nm x nm matrix is labeled by the vertices v of A* ^ which is helically wound on the torus 
(recall m = kn + 1). The vertex enumeration {1, 2, . . . , mn} can be made starting from the upper 
left corner following the southeast direction. Note that all vertices are visited in this way with the 
last site mn being neighbor to the first one. In addition, the four nearest neighbors of a vertex j 
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have the following coordinates (modulo mn): in the northeast and southwest directions, j"^^ = j + m 
and j*^™ = j — m, respectively; in the northwest and southeast neighbors, we have j"^™ = j — 1 and 

Here, the relevant N x N matrices are the identity I = In and the "forward shift" permutation. 



/ 1 ■ 
1- 



n = n 



N 



\ 



0- 

-10 0- 







(5.4) 



-n 



N-l 



which satisfies the following properties 11^ = — / and 11^ = 11 ^ 

We have chosen the shift matrix 11 with a phase e*'^ = — 1 in order to give odd parity to circuits 
winding around (see (HGj). 

From these definitions with N = mn, the matrix A can be written as 



A = /®f/ + n®x-n^®x^ + n™®r - (n^)" ® . 



(5.5) 



To block-diagonalize ()5.5p we introduce the N x N Fourier matrix Fn = F = {fij} defined as 
follows. Let uj = tUN = e"'^*/^, i = a/— 1, and note the following property 



We define 



fij '■- 



N 



(2i-l)j 



(5.6) 



(5.7) 



The Fourier matrix F is a unitary matrix, F F'^ = F"^ F = In, which diagonalizes the so-called 
class oi N X N "circulant" matrices. In particular, using ()5.7|) and property ()5.6|) . we have 

F-^UF = diag {uj, . . . , uj^^'^} = D (5.8) 

and F~^Il'^ F = D^. In addition, taking the Hermitian conjugate of ()5.8|) . gives 

[F-^UF)^ = F'^U^ F = D, 

where D = diag {u, u^, . . . , 

Now, we apply the Fourier matrix JF = Fmn ® h = F~^ ® J4) in ()5.5j) and use the 

multiplication rule ()5.1|) with ()5.8|) to obtain 

J^-M = / ® t/ + D ® X - D ® + ® F - L)™ ® = , (5.9) 



where 



Hi 



/ ua 

-ua 

\ Us -U3 - UJ^'"" 



—U4 — UJ^ —Us \ 

Uf U3 + LU^"^ 

-U5 

U5 J 



(5.10) 
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and j = 2(r — l)n + 2/ — 1 with r G {1, . . . , m} and / G {1, . . . , n}. 
As a consequence, the determinant of A can be written as 



detv4 = J]^ , (5.11) 



1=1 r=l 

and evaluated by the Laplace method, 

\Rj\ = — (a + 26cos (0^ + V'z') + 2ccosv9; + 2(icos (0^ + + 2ecos (0r — V'z")) ; (5-12) 
^1 

with 

(2r-l)7r (2/ - l)7r 

0r = , ^1 = , (5.13) 

m n 
ipi = Lfi/m + 7i/m, ipf = {1 ± m)(fi/m + Ti/m, and 

a = wl + W2+ wl + wj, 

b = Wi- 

c = W3-W2W4, (5.14) 

d = W3W4 — w^wq , 

e = W3W4 — W7 ws . 

Here, we have used = w^/wi and equation p.7|) to obtain this form. 

6 The Free Energy 

This section is concerned with the calculation of free energy / = /(/5), given by (jl.Hj) . of a class 
of quantum spin systems satisfying p.7p . Two limits are required to be taken: the Trotter limit 
m — *■ 00 of the partition function Zn^m and the thermodynamic limit n —>■ 00 of its logarithm. So 
far, we have obtained ^n,m exactly in Proposition 14.31 and equations ()5.11|) - ()5.14|) . for any value 
n, m G Z+. In the limit, these equations involve infinite products and the way we deal with them 
depends on whether m or n is taken to infinite. A third limit, lim^^oo /(/5); will also be taken in 
order to get the ground state energy. 

We begin by studying the Trotter limit. Taylor expanding ()3.2|) (with = 0) up to second 
order in 1/m and substituting the result into ()5.14|1 . gives 



c = —-^ + 



5 



(6.1) 



d = -1, 

^ = ~2 



where t] = — , n = + and ('^ = + rf. The value of d is exact in view of ()5.14j) and 
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We introduce a real valued function 



G{z) := -{a + c{z + z)) , 

and a complex valued function 

H{z) := 2^(1-"+™)/'" {d + bz + ez^) 
defined on the unit circle \z\'^ = 1 , and set 
A^{z) = G^-HH 

Ap^fh^ + e + K^ -2 ^2 



+ hK{z + z) + 



m" \ 2 4 

in view of (jHIl). Note that (G + A)(G - A) = HH. 

Using these definitions, equation ()5.12|) can be rewritten as 

- 2 Gl + Hi UJr + Hi UJr 



wl\Rj\ 



HifGi + Ai \fGi-Ai 



(6.2) 



(6.3) 



(6.4) 



(6.5) 



where Ur = e"^'', Gi = G'(e*'^') with (pr,fi as in ()5.13p and similar expressions holding for Hi and 
A,. 



Now, since {ur', r = 1, . . . , m} are the roots of —1, we have 

m 

Y[{X- OOr) = + 1 , 



(6.6) 



and the substitution of (|6.5|) into (|5.11|) yields 



+ 1 



(6.7) 



H {{-Hir + i-HiY'' + {Gi + AO" + (Gi - AO"} , 



1=1 



where we have used nr(~'^f) ~ 

To evaluate (j6.7j) . note that the asymptotic estimate 



m 



eM 1 + O 



m 



(6.8) 



holds for all c G M. From (j6.3p and the fact that m = kn + 1 (helical condition) for some positive 
integer k, we have 



[-Hi)"" = e-('=+^)"'^' (^l + O (J^^^ 
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which tends to 1 or —1 depending on whether k is odd or even. We always take k an odd number. 
From equations ()6.2|) and ()6.4p and asymptotic expansions (jHUI), we have 



iGi±Air = e^''^ il + Ol-] ] , (6.9) 



m 



where Si = 6{(fi) with 
This gives 

detA = n(2 + e-'+e--0(l + ^(^)) 



= n(2cosh<5,f + (6.11) 

= exp (^2 In (2 cosh 5,) j (l + o(l)) , 

and leads to the following expression for the partition function (see Proposition 14. 

Z„ = lim Z„,^ = exp |- / In [2 cosh (5(<^))]1 (l + o(-\\ , (6.12) 

where the limit is taken over the subsequence with m = kn + 1, k = 1,3,5,.... We have verified 
that the limit agrees with ()6.12|) if taken over the subsequence with m = kn, k G N+, corresponding 
to toroidal boundary condition. The remaining subsequences are expected to converge to the same 
limit. 

In approximating the Riemann sum by the Riemann integral in equation ()6.12|) we have used 
the fact that the integrand g = \ia{2 cosh 6) is a periodic and smooth function of (p to get the 
sharper than midpoint rule error estimate 



27r Jo 



for any G N+ (see Proposition 5.2 of jMFHj). 

From equations ()2.8|) and ()6.12|) . the free energy function ()1.3|) is thus given by 



^(^) = ^ 1^ [2cosh(5(v,))] (6.13) 

T^P Jo 



and the ground state energy, cq = lim /(/5), by 



eo = - [ dip ^h^ + (J^)2 + {jyy + 2h{.P - Jv) cosip + 2J^Jycos2ip , (6.14) 
Jo 

where we have used the fact that 5{ip) is symmetric about = tt. 

This concludes the analysis of the period one XY model in a transverse magnetic field. 
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7 XY Model with Period p = 2 



The previous analysis of a homogeneous quantum spin system with period p = 1 can be extended 
to arbitrary period. In this section we shall illustrate the procedure fully for the period p = 2 case. 

For a periodic model with period p, all steps of Section |21 can be repeated without modification. 
Since the purpose of that section was to separate consecutive terms of the Hamiltonian (jl.ll) we 
may adopt the same odd-even decomposition (j2.6j) . provided p is an even number. We also take 

= 2pn and impose periodic boundary condition. The only difference is that we now have p 
matrices p^^\ . . . , p^^^ , of the form ()2.14|) with distinct values (r^j, Kj, Q), i = 1, . . . ,p. 

Section ini is also maintained with a minor modification. A vertex v in the dual lattice A* ^ will 
be now distinguished according to its type labeled by {1, . . . ,p}- the label increases sequentially 
(modulo p) if one follows the northeast or southeast directions as seem in Figure 7. 



















V 




























(23 

V 



























Figure 7: Dual lattice with period p 

We associate to each vertex t>*^*) of type i a set {w^^\ . . . ,Wg^} of weights given by ()3.2|) with 
parameters {r]i,KiXi) and to each vertex (compatible) configuration ^ = {^v}, a weight w{^) = 
Yl^ w^^ with w^^ = w^^^ if V is of the type i. 

This completes the set up which allows Propositions 13.11 and 14.31 to be holden for any period p 
and we are now ready to evaluate the determinant of the correspondent matrix A. 

Up to this point there were no extra difficulties. The matrix A is defined by equation ()4.2|) with 
the signal determined by the orientation of the terminal graph Q = (A* as in Figure 6. But 
now the unit cells is composed by p different types of vertices which leads to 4p x 4p internal 
matrices. We shall write explicitly these matrices, whose rows and columns are indexed by the 
label of Figure 8, only for p = 2 . 

Similarly to what have been done in Sectional we set = ■ul^'' = w9'^ /w^^ and u'^ = = 
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Figure 8: Unit cell for period p = 2 



w'p /wf\ and let 



V 



and 










—U4 


-Ms 














\ 




-Me 







U3 




















-M7 





-Ms 


1 















Ms 


-M3 


Ms 



























-1 








u'q 


-M4 - 


-m's 























u'j 




















M4 


—u'j 


- 






v 














m's 




< 





/ 


Q 


-{ 


^ 
X , 


)■ 


R = 


\ 


j 


and 




s 



u 
-x^ 



X 
U' 



(7.1) 




Y 



(7.2) 



where U, X and y is as in ()5.2j) and ()5.3p . If m = /cn + 1, with an odd number, we recall that 
A* ^ becomes helically wound on a torus. The matrix A can thus be written as (see Figure 9 for a 
self-explanation) 



A = I ®v + n®Q - If + ®R- [n^)"^ ® + ® s - [if)'^^® , (7.3) 

where / and 11 are the nm x nm identity and the forward shift. 




Figure 9: Numbering nearest neighbors unit cells 

Conjugating A by the Fourier matrix JF = Fmn ® h gives 

U 



-Aft 
U' ' ' 



(7.4) 
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where j = 2(r — l)n + 21 — 1, with r G {1, . . . , m}, / G {1, . . . , n} and 
with uj = e™l . This imphes 











-1 





\ 













^i(m+l) 



















v 













/ 



(7.5) 



where 



detA=nni^: 

l=\ r=l 

15-^.1 = |f/| U' + MjU~^M^. 



(7.6) 
(7.7) 



Here, we used the fact that ()7.4p can be brought to block diagonal form by a non-unitary trans- 
formation 



1 

-MjU-^ 1 



U -M. 
M, U' 



1 u-^m\ 



U 

u' + m,u-^m] 



Using Laplace method and equation ()4.inj) . we have \U\ = u\ and 





f 





-Ms 


-Ms 


M7 
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1 




Ms 





-Ms 


M4 








Ms 


Ms 





Me 






\ 


-My 


— M4 


-Me 





/ 



which, in view of fl7.5p and ()7.7|) . gives 



— —9 

Mi 



U2U'^ - ^6^"^'^"'+^^ 



— M2M4 — MsCJ 



M2M4 + UsCo'-'' 





-M2M7 + UyCiJ 



M2M7 - MyCU^'™ M2M3 + M4C^^'(2'™+^) 







Using Laplace method and ()4.10|) . equation ()7.8|) can be written as 
{wiw[f\Sj\ = 2a + 26cos(^j) + 2ccos(m^j) + 2rfcos((m + 1)^^) 

+2ecos(2m^j) + 2/cos((2m + 1)^^) + 2^cos(2(m + 1)^^) 
where = irj/nm = (pr + fi/Tn — it /m. Defining 

Ql = WiW[ + ^2^2 5 
^2 = WsW'^ + WiW'^ , 
fig = W^w'q + WqW'^ , 



(7.8) 



(7.9) 



(7.10) 
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the a-g constants are given by 

2a = Ql + Ql + Ql + Ql + 2{w5WQWjWg + w-jW^w'^Wq — wiW2w[w'2 — wsw^w'^w'^) , 

d = fi2^^4 - ^^lf^3 , (7.11) 

e = - WyWs) (tf 3W4 - ^7^) ' 

9 = (W3W4 - w^we) (ti'>4 - «) , 

Note that, in view of equations in ()3.2p . b = f and g = 1 hold exactly. This fact together with 
the trigonometric relations 

cos(^j) + cos(2m^j + Q = 2 cos((m + 1)^^) cos(m^j) , 

and 

cos (2(m + 1)Q = 2 cos^((m + 1)^^) - 1 , 
lead the right hand side of (j7.8p to be written as 

{wiw[y\Sj\ = A {cos^ {{m + 1)Q -2Bicos{{m + l)^j) + Ci) (7.12) 

= 4 (cos((m + - xl) (cos((m + 1)^^) - x^) , 

where 

Bi = -l- + -cos(^,, 

^ a - 1 c / X e , , 
Ci = + - cos(v?/) + - cos{2ipi) 

(note that cos(m^j) = cos(2(r + l)7r + (2/ — l)/n) = cos{ipi)) and 



x^ 



Bi±^Bf-Ci (7.13) 
= l + ^{/€2 + (/€')' + ^' + (^')' + (/^ + /^')' + 2(/«/s:'-W)cos(^0}± A + , 

where 

Df = ^{ [f^V - + i^^v' - ^^'v) cosiifi)]^ + ih + h'f [k" + {k'Y + 2kk' co^i^i)] } . (7.14) 

Here a word about the 1/m expansion is in order. Differently from the period p = 1 case, the 
constants ()7.11|) of period p = 2 have to be computed up to order 1/m'' to pick up the relevant 
terms. This happens because B^ — C, whose leading order term is written in equation ()7.14|) . is 
O (1/m^). So the square root Di has the same order 1/m^ of 5/ — 1. Equation (j7.12|) reduces the 
present problem of evaluating det A to the one already dealt with in the previous section. 
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Repeating the steps in ()6.2|) - ()6.5p . we have 



(7.15) 

where Hi = (a;;)-(i-'^+™)/™ and (Af )^ = (x^)^ - = (xf )2 - 1 which, in view of ^J^, can 
be estimated by 

(A±)' = 4 l'^' + ('^')' + + Wf + ih + h'f + (/t/t' - r,i) cos{vi)] ±2Di + (—\ . (7.16) 

Note that Af = C (1/m). 

Now, repeating the steps from (jfj.fij) to ()fi.ll|l . leads ()7.fjj) to be written as 

(wiw;i)MetA = exp [ln(2cosh(5,+) +ln(2cosh5f)] j + O (^^^ , (7.17) 

where 5f = ^^(v^/) = hm mAf/2 which, when expressed in terms of the couphngs and magnetic 

m— >oo 

fields, is given by 

{6H^)f = ^ {{Jtf + {r,f + {Jlf + {jyf + 2{JIJI + Jf J^^') cosy.) 

±Y { [(^n' - {J2f - {Jlf + (^1)' + 2(Jf J| - Jf J^^') cosy.]' 

Hhi + /i2)' [(Jf + Jlf + (^2 + ^D' + 2(Jf + Jf J| + JIJ2 + </i ^1) cos<^] . 
The free energy function (jl.Hj) for p = 2 periodic systems is thus given by 

/2(/?) = r(/5) + r(/?), (7.18) 

with 

f^iP) = l^[d^ In [2cosh(5±(y.))] . (7.19) 
Note that the total size of the system is now equal to 4n and this explain why ()7.19|) is formally 

1/2 of dnm. 

It is interesting to observe that, if we set Jf = Jf = J^, Jl = Jl = J^ ^"^^ hi = h2 = h, we 
have 

{6^{ip)f = /3 {{rf + {.pf + 2.rjy cosy. ± 2h{r + .p) |cos((/./2)|) , 

and this gives (since 5^(y) is symmetric about y? = tt) 



/2(/5) 



I / j-K/2 j-TT/2 \ 

— / rfrln [2cosh(5+(2T))] + / dr In [2 cosh (r (2r))] , (7.20) 



whose combination of the two terms yields (|6.13|) . 
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